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ABSTRACT 

We present Herschel PACS photometry of seventeen B- to M-type stars in the 30 Myr-old Tucana- 
Horologium Association. This work is part of the Herschel Open Time Key Programme "Gas in 
Protoplanetary Systems" (GASPS). Six of the seventeen targets were found to have infrared excesses 
significantly greater than the expected stellar IR fluxes, including a previously unknown disk around 
HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate 
the temperatures and abundances of the dust in the systems. For the five stars that show excess 
emission in the Hershcel PACS photometry and also have Spitzer IRS spectra, we fit the data with 
models of optically thin debris disks with realistic grain properties in order to better estimate the disk 
parameters. The model is determined by a set of six parameters: surface density index, grain size 
distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The 
best fitting parameters give us constraints on the geometry of the dust in these systems, as well as 
lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally 
resolved PACS 70 /zm imaging. 

Subject headings: stars: circumstellar matter — infrared: stars 



1. INTRODUCTION 

Debris disks are the last stage of circumstellar disk 
evolution, in which the gas from the protoplanetary and 
transitional disk phases has been dissipated and the dust 
seen comes from collisions between planetesimals. In the 
youngest debris disks, < 100 Myrs old, terrestrial planets 
may still be forming (Kenyon & Bromley 2006). Giant 
planets must form before the gas dissipates; their grav- 
itational interactions with planetesimals and dust can 
leave signatures in debris disks. Cold debris disks may 
be Kuiper-belt analogs, signaling the location and prop- 
erties of planetesimals remaining in the disk. 

To be suitable for life, terrestrial planets in habit- 
able zones must have volatiles such as water brought to 
their surfaces from beyond the ice line. Planetesimals in 
Kuiper Belt-like debris disks may provide this reservoir 
of volatiles (Lebreton et al. 2011). We do not yet have 

1 Department of Astronomy, University of Maryland, College 
Park, MD 20742; jessd@astro.umd.edu 

2 Exoplanets and Stellar Astrophysics Laboratory, NASA 
Goddard Space Flight Center, Code 667, Greenbelt, MD 20771 

3 Space Telescope Science Institute, 3700 San Martin Dr., Bal- 
timore, MD 21218 

4 UJF - Grenoble 1 / CNRS-INSU, Institut de Planetologie 
et d'Astrophysique de Grenoble (IPAG) UMR 5274, Grenoble, 
F-38041, France 

5 ALMA, Avda Apoquindo 3846, Piso 19, Edificio Alsacia, Las 
Condes, Santiago, Chile 

6 European Southern Observatory, Alonso de Cordova 3107, 
Vitacura, Santiago, Chile 

7 Dpt. Fi'sica Teorica, Facultad de Ciencias, Universidad 
Autonoma de Madrid, Cantoblanco, 28049 Madrid, Spain 

8 Astrophysikalishes Institut, Friedrich-Schiller-Universitat 
Jena, SchillergaBchen 2-3, 07745 Jena, Germany 

9 Institute for Astronomy (If A), University of Hawaii, 2680 
Woodlawn Dr., Honolulu, HI 96822 

10 Centra de Astrobiologi'a Depto. Astrofisica (CSIC-INTA), 
POB 78, 28691 Villanueva de la Canada, Spain 

11 SOFIA-USRA, NASA Ames Research Center, Building 
N232, Rm. 146, Moffett Field, CA 94035 



the capability to detect the planetesimals in these disks, 
but we can detect the smaller dust grains. These dust 
grains are believed to be produced through the collisions 
of the larger planetesimals, and therefore are likely to 
have similar compositions to the larger undetected bod- 
ies. The properties and locations of the dust grains in 
Kuiper Belt analogs can provide clues to the properties of 
the hidden planetesimals. Additionally, dust structures 
in the disk may point to unseen exoplanets (e.g. /3 Pic b; 
Lagrange et al. 2010). 

To measure the cold dust in the outer regions of de- 
bris disks, we need great sensitivity at far-infrared wave- 
lengths where the thermal emission from the cold dust 
grains peaks (> 70 /im). The Herschel Space Observa- 
tory (Pilbratt et al. 2010) provides a unique opportunity 
for sensitive debris disk surveys. HerscheFs PACS in- 
strument (Poglitsch et al. 2010) is sensitive to the cold 
dust with a wavelength range of 55 — 210 ^tm. Addition- 
ally, HerscheFs spatial resolution is almost 4 times better 
than Spitzer at similar wavelengths, and therefore reduc- 
ing confusion with background galaxies and interstellar 
cirrus and making it easier for Herschel to detect faint, 
cold debris disks. 

In this paper, we present results of a sensitive Herschel 
debris disk survey in the 30-Myr-old Tucana-Horologium 
Association. This work is part of the Herschel Open 
Time Key Programme "Gas in Protoplanetary Systems" 
(GASPS; Dent et al. in prep, Mathews et al. 2010). The 
GASPS survey targets young, nearby star clusters with 
well determined ages, ranging from 1-30 Myrs. This 
range of ages covers the stages of planet formation from 
giant planet formation (~ 1 Myr; Alibert et al. 2005) 
to the late stages of terrestrial planet formation (10-100 
Myrs; Kenyon & Bromley 2006). The targets in each 
group were chosen to span a range of spectral types from 
B to M. 

The 30 Myr-old Tucana-Horologium Association is the 
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TABLE 1 

Stellar Properties in the 30 Myr-old Tucana-Horologium Association 
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a Spectral types listed are from the SIMBAD Astronomical Database 
b Calculated from stellar modeling. See Section 3 

c Distances are taken from the Hipparcos Catalog (Perryman & ESA 
d Binary or multiple star system 
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oldest in the GASPS survey. The Tucana-Horologium 
Association, discovered independently by Zuckerman & 
Webb (2000) and Torres et al. (2000), is a group of ~ 
60 stars with common proper motion and an average 
distance of 46 pc (Zuckerman & Song 2004). About ~ 
1/3 of the targets have debris disk systems known from 
previous Spitzer surveys (Hillenbrand et al. 2008; Smith 
et al. 2006). 

We obtained Herschel PACS photometry of the sev- 
enteen GASPS targets in the Tucana-Horologium Asso- 
ciation. We also obtained PACS spectra for two of the 
targets. Previously unpublished Spitzer IRS spectra for 
three targets are presented. In Section 2, we present our 
methods and results of data reduction and aperture pho- 
tometry. In Section 3, we fit blackbody and modified 
blackbody models to the detections and upper limits to 
determine dust temperatures and fractional luminosities. 
We further analyze some of the disks with our optically 
thin dust disk model in Section 4. Additionally, we dis- 
cuss the detection of a marginally resolved disk in our 
sample in Section 5 and present conclusions in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

The stellar properties of the seventeen Tucana- 
Horologium stars observed are listed in Table 1. These 
stars have spectral types ranging from B9-M3 and dis- 
tances of 25-64 pc. PACS scan map observations were 
obtained for all seventeen Tucana-Horologium targets at 
70 and 160 fim. Additional observations at 100 /im were 
taken for five targets. The stars were observed at two 
scan angles, 70deg and 110 deg. The two scans were 
combined to reduce the excess noise caused by streaking 
in the scan direction, as suggested by the Herschel PACS 
Instrument Calibration Centre (ICG) 12 . Each scan map 
was executed with the medium scan speed (20"s ) and 
consisted of 10 legs with scan lengths of 3' and scan leg 
separation of 4". The number of repetitions varied from 
1 to 4 based on the expected flux density of the star. The 
on-source time for each observation is given in Table 2. 



Two targets, HD105 and HD3003, were also observed in 
the PACS LineScan mode at 63 fim and 190 /im, target- 
ing the OI fine structure line and DCO+ respectively. 

2.1. Photometry 

The data were reduced with HIPE 7.2 (Ott 2010) us- 
ing the standard reduction pipeline. The final maps have 
pixel scales of 3.2/3.2/6.4"pixel _1 in the 70/100/160 /im 
images respectively, corresponding to the native pixel 
scales of the PACS detectors. The two scans were re- 
duced separately, then averaged together. The flux 
values were measured using an IDL aperture pho- 
tometry code with apertures of 5.5/5.6/10.5" for the 
70/100/160 /im images respectively, as recommended for 
faint sources. Aperture corrections were applied based 
on the encircled energy fraction from PACS observations 
of Vesta provided by the Herschel PACS ICC 13 . The sky 
annulus for error estimation was placed 20-30" away from 
star center for the 70 and 100 /im images and 30-40" away 
for the 160 /im images. For three of the targets, HD2884, 
HD2885 and HD53842, the fields were contaminated by 
nearby stars and background galaxies, so the sky annulus 
was offset to a nearby clean field. 

The RMS pixel uncertainty, a P i X , was estimated by 
calculating the standard deviation of the pixels in the sky 
annulus. The total statistical error in the measurement 
is given by 



Ctot 




(1) 



* corr~ lj corr 



where a corr is the aperture correction, x corr is the cor- 
related noise correction factor 13 (0.95 at 70 and 100 /im, 
0.88 at 160 /im), and n ap and n s k y are the number of 
pixels in the aperture and annulus respectively. An ab- 
solute calibration error was also added in quadrature 

12 PICC-ME-TN-036: http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/PhotMiniScan_ReleaseNote_20101112.pdf 
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TABLE 2 

Herschel PACS Photometry Results 





70 /im 


100 /urn 


160 fim 


Star 


Obs. Flux* 
(mjy) 


Stellar Flux 

(mjy) 


On-source 
Time (s) 


Obs. Flux 
(mjy) 


Stellar Flux 
(mjy) 


On-source 
Time (s) 


Obs. Flux 
(mjy) 


Stellar Flux 
(mjy) 


On-source 
Time (s) 


Debris Disks 


HD105 


128.3 ± 7.0 


3.1 


72 


151.2 ±5.7 


1.7 


144 


81.2 ± 12.3 


0.6 


144 
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59.7 ±3.8 
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72 


19.0 ±2.4 


4.1 
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13.0 ±0.9 
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0.6 
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3.2 
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0.6 
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360 
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0.3 


360 
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10.5 ±0.8 


2.3 
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< 15.5 


0.5 


720 


Non-Excess Stars 


HD2884 


7.7 ± 1.4 


12.1 


360 








< 18.6 


2.3 


360 


HD 16978 


15.4 ± 1.7 


15.8 


144 








< 14.4 


3.1 


144 
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8.6 ±1.1 


8.7 


540 








< 10.0 


1.7 


540 


HD2885 


16.0 ± 1.4 


16.5 
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360 


HD53842 


9.9 ±3.3 


2.3 
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< 9.2 


1.3 


144 
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0.5 


144 


HD44627 


< 3.5 
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360 
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0.3 


360 
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360 
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360 


HD3221 


< 3.7 


2.4 


360 
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0.4 


360 


HIP107345 


< 3.3 


0.9 


360 








< 10.8 


0.2 


360 


HIP3556 


< 16.5 


1.3 


72 


< 16.9 


0.7 


144 


< 17.4 


0.2 


144 


GSC8056-482 


< 3.6 


1.6 


360 








< 20.4 


0.3 


360 



: la uncertainties include statistical and 
were calculated by fitting the optical and 



ystcmatic errors added in quadrature. The expected stellar fluxes 
near-IR data with NextGen stellar atmosphere models (Hauschildt 



at the PACS 
ct al. 1999). 



wavelengths 



with the statistical error to give the total uncertainty re- 
ported in Table 2. The absolute calibration errors given 
by the Herschel PACS ICC are 2.64/2.75/4.15% for the 
70/100/160 /im images respectively 13 . Upper limits for 
all non-detections were determined using the 3<7 errors 
from the same aperture photometry method as above. 

The PACS photometry fluxes and upper limits are 
listed in Table 2. Six of the seventeen targets were deter- 
mined to have infrared excesses above the photosphere, 
including one previously unknown disk, HD30051. A sig- 
nificant excess is defined here as a 70 /im photosphere- 
subtracted flux greater than 3a. One of our targets, 
HD2884, was previously suspected to have an infrared 
excess (Smith et al. 2006). In the PACS data, this tar- 
get seems to have less flux than the expected photo- 
spheric value at 70 /im. The Spitzer IRS spectrum also 
shows no hint of an IR excess. Smith et al. (2006) notes 
the possibility of contamination of HD2884 with a back- 
ground galaxy, but could not separate the two sources 
with Spitzer's spatial resolution. Figure 1 shows the con- 
fused field of HD2884 observed with Herschel and Spitzer 
at 70 /im. HerscheFs improved spatial resolution allows 
us to properly resolve the sources and avoid contami- 
nation. We determine that HD2884 has no detectable 
debris disk. 

HD53842 has a photosphere subtracted flux that is 
only 2.3ct above the noise at 70 [im. The Spitzer IRS 
spectrum from Moor et al. (2009) also shows a marginal 
excess at 14-35 /im but may suffer from contamination. 
The Spitzer 24 /im image shows two nearby stars within 
15" whose Point Spread Functions (PSFs) appear to 
overlap with that of HD53842. These sources are also 
seen in the PACS images (see Fig. 2) and make an accu- 
rate determination of the flux difficult. We are unable to 
confirm if HD53842 has a debris disk. 



13 PICC-ME-TN-037: http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/pacs_bolo_fluxcal_report_vl.pdf 



2.2. Spectroscopy 

Two targets, HD105 and HD3003, were observed in the 
LineScan mode of the PACS instrument. The lines tar- 
geted were [OI] at 63.185 /mi and DCO+ at 189.570 /mi. 
The continuum detection limit was 202 mjy at 63 /jm. 
The photometry detection at 70 /im was below this limit 
for both HD105 and HD3003, and the noise levels in 
the spectra are comparable to or larger than the flux at 
70 /xm. Therefore, we conclude that the continuum is 
not detected in these observations. There are also no 
emission lines detected above the noise. 

Upper limits to the line fluxes were calculated by in- 
tegrating over a Gaussian with a width equal to the 
instrumental Full- Width Half-Maximum (FWHM) and 
a height given by the standard deviation of the noise, 
placed at the wavelength where the line is expected. The 
upper limits are given in Table 3. 

The three targets with previously unpublished Spitzer 
spectra (HD1466, HD2884 and HD3003) were observed 
with the IRS spectrograph (Houck et al. 2004), using 
the Short-Low (5.2-14 /im) and Long-Low modules (14- 
38 /im; A/AA ~ 60). These three systems were observed 
as part of the Spitzer GTO program 40651 (PI: J. Houck). 
Since the sources were expected to be unresolved and 
the pointing of Spitzer was excellent, all of our targets 
were observed in Staring mode with no peak-up. We 
carried out the bulk of the reduction and analysis of our 
spectra with the IRS team's SMART program (Higdon 
et al. 2004). 

3. BLACKBODY AND MODIFIED BLACKBODY FITS 

We constructed spectral energy distributions for all 
sources using data from Hipparcos, 2MASS, Spitzer 
(IRAC, IRS, and MIPS), AKARI, IRAS, the WISE pre- 
liminary release, and our new Herschel data. The SEDs 
of debris disk stars are shown in Figure 3 and those of 
non-excess stars in Figure 4. The data used for these 
figures are listed in Table ?? (online version). For each 
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Fig. 1 — The field of HD2884 and HD2885 at 70 /im with Her- 
schel (top) and Spitzer (bottom). The better spatial resolution 
of Herschel allows the sources to be cleanly separated and avoids 
contamination of photometry. 



TABLE 3 

Herschel Spectroscopy: 3<t Line Upper Limits 



Star 


Line 


Wavelength 


3<r Upper Limit 






(/urn) 


(W/m 2 ) 


HD105 


[OI] 


63.185 


< 9.63 x 10~ 18 




DCO+ 


189.570 


< 1.18 x 10~ 17 


HD3003 


[OI] 


63.185 


< 1.43 x 10~ 17 




DCO+ 


189.570 


< 8.55 x 10~ 18 



target, we fitted the stellar photosphere with NextGen 
models of stellar atmospheres (Hauschildt et al. 1999). 
The best fitting stellar model was determined through 
X 2 minimization with effective temperature and the nor- 
malization factor as free parameters. The Spitzer IRS 
spectra were binned to a resolution of AA/A ~ 0.1 and 
only the data with A > 8 pim were used to determine the 
fit of the excess. 



Fig. 2.— The field of HD53842 at 70 fim with Herschel. The star 
is barely detected at 70 ^m (3cr detection). The nearby sources 
may have contaminated the Spitzer IRS spectrum of HD53842. 

For our six disk detections, we fitted the excess emis- 
sion with a standard single temperature blackbody model 
with two parameters, the temperature and the fractional 
disk luminosity. This model is a simplified disk represen- 
tation that assumes all the dust is at the same temper- 
ature, and behaves like a perfect blackbody. While this 
model is simple, it gives a good first estimate of the dust 
temperature and abundance, and unlike more physical 
models, one can get a fit even with very few data points. 
The best fit was determined through \ 2 minimization. 
We also calculated the goodness-of-fit for each best fit- 
ting model for comparison with more detailed models 
described in Section 4. The goodness-of-fit, Q, is the 
probability that a x 2 value this poor will occur by chance 
given the error in the data. Q is given by the normalized 
incomplete gamma function 

\ T(a,x) N x 2 

where T(a,x) and T(a) are incomplete and complete 
gamma functions respectively, and N is the number of 
degrees of freedom. Q varies from to 1, with larger 
values indicating better fits. An acceptable model is one 
that has a value of Q > 1(T 3 (Press et al. 1992). 

The Q values, however, are largely influenced by the 
IRS spectrum. The IRS spectra have more data points 
compared to the few PACS photometry points. There- 
fore, fits that may appear to be just as good, can have 
very different Q values if the uncertainties in the IRS 
spectra are different. The point of the Q value is not to 
compare the goodness-of-fit between disks, but to com- 
pare between models of a given disk. 

From the fits, we determined the disk temperatures 
and the fractional infrared luminosities (Lir/L*). For 
the remaining sources, we determined upper limits on 
Lir/L* by fitting blackbody models to the flux upper 
limits. Upper limits on Lir/L* depend on temperature, 
so the largest Ltr/L* value found assuming dust tem- 
peratures from 10-300 K was adopted. The results of the 
fits to the disks are given in Table 4. The disks display a 
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Fig. 3. — Single temperature blackbody fits to the SEDs of the 6 debris disks in the sample. Filled circles are data from the literature, 
open squares are the Herschel PACS data. The solid black curve is the fit to the stellar photosphere using the NextGen models (Hauschildt 
et al. 1999). The red curve shows the best fitting single temperature blackbody model. The figure label for each disk displays the best fitting 
temperature and the fractional IR luminosity Lip^/L*. The dashed red lines show the best-fitting modified blackbody models. Available 
MIPS 70 fj,m data are not plotted here, as they are consistent with the PACS data, which have smaller uncertainties. 



large range of temperatures and Lir /D* values, showing 
no correlation with spectral type. 

We also fit the disks with a modified blackbody model. 
This model assumes that at longer wavelengths, the 
grains no longer emit as blackbodies, but instead the dust 
emissivity is decreased by a factor (A/A ) _/3 for A > Ao, 
due to a lower dust opacity at longer wavelengths com- 
pared to a blackbody. Here Ao = 2irabiow, where auow 
is the blowout size due to radiation pressure. We do not 
have sub-mm or mm data to constrain the parameter /3, 
therefore we assume j3 = 1.0 for these models based on 
previous measurements for debris disks (Dent et al. 2000; 
Williams & Andrews 2006; Nilsson et al. 2010). ISM dust 
has (3 k, 2.0 (Boulanger et al. 1996), but debris disks are 
expected to have lower values of /3 due to the presence 
of larger grains (Draine 2006). The results of the modi- 
fied blackbody fits are also shown in Table 4. The use of 
modified blackbody models shifts the characteristic tem- 
peratures to lower values, a trend which was also noticed 
by Carpenter et al. (2009). The goodness-of-fit was im- 
proved for HD105, HD3003, and HD12039 with the use 
of the modified blackbody model, but made worse for 
HD1466 and HD202917. 

The blackbody and modified blackbody models pro- 
vide a first look at the disk properties. For example, 
disks with low characteristic temperature such as HD105, 
HD30051, and HD202917, are likely to have large inner 
gaps and copious amounts of cold dust far from the star. 
To test this idea, we need a more physical model of the 
disk. 



4. DUST DISK MODEL 

To further investigate the disk properties, we fit the 
disks with an optically thin dust disk model. Rather 
than assuming the grains to be perfect blackbodies, the 
dust model assumes the grains have a particular emissiv- 
ity that is dependent on the size of the dust grains and 
the wavelength of radiation being absorbed or emitted. 
For the Tucana-Horologium disks, we make the assump- 
tion that the grains are purely silicate in composition, 
specifically, astronomical silicates (Draine & Lee 1984; 
Laor & Draine 1993; Weingartner & Draine 2001). The 
poorly populated SEDs and lack of resolved imaging pre- 
vents us from fitting more complex grain compositions. 

The geometry of the disks is described by three param- 
eters: the inner and outer radius (r m in and r max ) and 
the radial surface density profile. The azimuthally sym- 
metric radial surface density profile is characterized by 
a power-law with index q, such that the surface density 
E varies with radius as S(r)dr ex r q dr. A power-law in- 
dex of q — is expected for a transport dominated disk, 
and an index of q = —1.5 is expected for a collisionally 
dominated disk (Krivov et al. 2006; Strubbe & Chiang 
2006). 

The population of dust grains is also characterized 
by three parameters: the minimum and maximum dust 
grain sizes in the disk (a m i n and a ma x), and the dis- 
tribution of grains with sizes between a m i n and a max . 
The grain size distribution is typically assumed to be a 
power-law with index k such that the number density 
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Fig. 4. — Spectral energy distributions for the 11 stars that were not found to have debris disks. The solid curve is the fit to the stellar 
photosphere using the NextGen models (Hauschildt ct al. 1999). Upper limits for Ljr/L* are shown in Table 4. Available MIPS 70/j.m 
data are not plotted here, as they are consistent with PACS data for every source except HD2884, which lies in a confused field (see Fig. 
1). HD53842 has only a marginal detection of an excess at 70 fim and also lies in a confused field (see Fig. 2). Consequently, wc arc unable 
to determine if HD53842 has a debris disk. 
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TABLE 4 
Model Fit Results 



Debris Disks 



Blackbody Fits 



Star 


Temp (K) 


WL* 






Q 


HD105 


49.7 ±1.4 


(2.59 ±0.05) x 


10" 


-4 


0.004 


HD3003 


207.9 ±2.5 


(1.21 ±0.01) x 


10" 


-4 


0.00 


HD1466 


145.5 ±6.7 


(7.21 ±0.11) x 


10" 


-5 


0.61 


HD30051 


57.8 ±6.4 


(3.27 ±0.20) x 


10" 


-5 


0.82 


HD202917 


80.5 ± 1.5 


(2.50 ±0.01) x 


10" 


-4 


0.93 


HD12039 


111.6 ±6.2 


(6.55 ±0.39) x 


10" 


-5 


0.38 


Modified Blackbody Fits 


HD105 


40.3 ± 1.2 


(2.41 ±0.90) x 


10" 


-4 


0.11 


HD3003 


160.6 ± 1.6 


(1.08 ±0.01) x 


10" 


-4 


0.00 


HD1466 


106.8 ±4.3 


(6.33 ± 0.03) x 


10" 


-5 


0.002 


HD30051 


39.8 ± 2.8 


(2.81 ±0.05) x 


10" 


-5 


0.82 


HD202917 


64.6 ± 1.1 


(2.50 ±0.01) x 


10" 


-4 


0.13 


HD12039 


85.1 ±3.7 


(6.33 ±0.33) x 


10" 


-5 


0.43 


Non-Excess Stars 




Blackbody 


Modified Blackbody 






Lir,/L„ 


Lir/L* 








HD2884 


< 2.9 x 10~ 6 


< 6.7 x 10" 


-6 






HD16978 


< 1.9 x 10~ 6 


< 8.6 x 10" 


-7 






HD224392 


< 5.0 x 10~ 6 


< 1.6 x 10" 


•5 






HD2885 


< 2.7 x 10~ 6 


< 1.8 x 10" 


-5 






HD53842 


< 6.4 x 10~ 5 


< 2.9 x 10" 


-4 






HD44627 


< 1.9 x 10~ 4 


< 1.7 x 10" 


-4 






HD55279 


< 8.6 x 10~ 4 


< 2.5 x 10" 


-3 






HD3221 


< 2.4 x 10~ 4 


< 1.3 x 10" 


-3 






HIP107345 


< 8.4 x 10~ 4 


< 4.2 x 10" 


-4 






HIP3556 


< 1.2 x 10~ 3 


< 5.4 x 10" 


-4 






GSC8056-482 


< 1.2 x 10" 3 


< 5.3 x 10" 


-4 







varies with grain size a as n(a)da oc a K da, where a value 
of k ~ —3.5 is expected for a steady-state collisional cas- 
cade (Dohnanyi 1969). 

The model iteratively determines the equilibrium dust 
temperature at each radius for each grain size by bal- 
ancing radiation absorbed from the star with radiation 
re-emitted through the formula 

f°° / R \ 2 

J Qabs{v,a)\-y \ B v (T eff )dv = 

/ 4Q abs (is,a)B v (T d (r,a))dv, (3) 
Jo 

where B u (T e ff) is the blackbody flux coming from the 
surface of the star with temperature T e f f and radius i?* , 
B u {Td{r, a)) is the blackbody flux radiating from the dust 
grain of size a at a distance r away from the star, and 
Qabs{ v i a ) is the dust absorption coefficient calculated by 
Draine & Lee (1984) for an astronomical silicate grain of 
size a. The total flux is then determined by summing up 
over all radii and grain sizes according to the formula 

F„=A f""" (—Yda 

Jo m ,„ \ a max / 

■[ 47ra 2 Q Q6s (^a)7rB ly (T d (r,a))27rr (-^— ) dr, 

(4) 

where A is a normalization constant that includes the 
distance to the system and the total amount of material 
in the disk. 



4.1. Model Parameters 

The Tucana-Horologium debris disks are faint (all are 
under 150 mJy at 70 /an), so the SEDs of these targets 
are not well sampled due to lack of detections at the 
longer wavelengths. For this reason, we must fix some 
of the model parameters. In particular, because of the 
lack of sub-mm and longer wavelength data, we are un- 
able to constrain parameters that affect this region of 
the SED, specifically cold grains that are larger than 1 
mm or farther out than ~ 120AU. Therefore, we fix the 
maximum grain size and the outer radius to these val- 
ues. For HD202917, we fixed the outer radius to 80 AU 
to be consistent with HST imaging (Krist 2007; Mustill 
& Wyatt 2009). The radial surface density profile also 
cannot be constrained with the current data. Hence we 
fix the power-law index of the radial density profile to 
q = —1.5, the value expected for collisionally dominated 
disks. 

The SEDs of HD105, HD3003, and HD202917 are pop- 
ulated enough to get well constrained fits with only these 
three parameters fixed. However, HD1466 and HD12039 
were only detected at 70 /im with PACS. The models for 
these disks were further constrained by fixing the min- 
imum grain size and the grain size distribution. The 
grain size distribution was fixed to a power-law index of 
k = —3.5, the value expected for a steady-state colli- 
sional cascade, and the minimum grain size was fixed to 
the expected blowout size for astronomical silicates. This 
is the size at which the radiation pressure is half of the 
gravity; grains this size and smaller are ejected from the 
system. The blowout size is calculated from the formula 
given in Backman & Paresce (1993) as 

assuming constant density, spherical grains produced 
from planetesimals on circular orbits. This equation as- 
sumes that the grains are of a certain composition and 
all have the same density, and that the density is uniform 
throughout the grain. For astronomical silicates, we as- 
sume a grain density of 2.5 g cm -3 . This simple equation 
provides a good estimate of the expected minimum grain 
size, but it relies heavily upon many assumptions and 
must be treated solely as an estimate, accurate perhaps 
to only an order of magnitude. 

4.2. Results 

To determine the best model for each system, we com- 
pared a grid of models to the SED and found the set of 
parameters that produced the x 2 minimum value. The 
error of each fit parameter was determined from the la 
confidence interval in the \ 2 distribution after fixing the 
other parameters to their best fit values. The best fitting 
parameters and the errors are displayed in Table 5, and 
the best fitting SED models are shown in Figure 5. The 
HD30051 disk could not be fit at all with this model be- 
cause its SED has too few data points, having not been 
known to have a debris disk before Herschel. 

Lower limits on the disk masses are also given in Table 
5. These masses are calculated by summing up the mass 
of silicate grains with an assumed uniform density of 2.5 
g cm" 3 , with sizes between a m j n and a max and disk radii 
between r m ; n and r max . The mass also depends on the 
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Fig. 5. — Best fitting models for the five of the debris disks with well populated SEDs. The solid black lines shows the stellar photospheres 
that are fit with the NextGen models (Hauschildt et al. 1999). The best fitting models of the disks are shown in red. The best fitting 
model parameters are given in Table 5. 



normalization constant and the radial density profile of 
the disk. It is a lower limit on the mass because it only 
takes into account the dust smaller than 1 mm, and not 
larger pebbles and planetesimals hidden in the disk. 

4.2.1. ED 105 

HD105, a GOV star 40 pc away, has the largest 
IR excess in the PACS wavebands of all the Tucana- 
Horologium disks in this sample, yet shows very little 
mid-IR excess in the Spitzer IRS spectrum. This disk 
was observed with Spitzer as part of the Formation and 
Evolution of Planetary Systems (FEPS) Legacy Survey 
(Meyer et al. 2006). Meyer et al. (2004) first fit the disk 
with models from Wolf & Hillenbrand (2003) that used 
Draine & Lee (1984) grain emissivities for astronomical 
silicates and graphite in ISM ratios and found an inner 
radius of 45 AU for a fixed minimum grain size of 5 /an. 
They assumed a flat radial density profile (q=0), a grain 
size distribution power-law index of k = —3.5, an outer 
radius of 300 AU, and a maximum grain size of 100 /im. 

Hollenbach et al. (2005) found an even lower x 2 value 
using the Wolf & Hillenbrand (2003) models with a min- 
imum grain size of 21 /im and a lower inner radius of 19 
AU, but confirmed that the inner cavity must be > 13 
AU because of the lack of IR excess at A < 35 /im. The 
Spitzer IRS spectrum also shows no evidence of gas lines, 
and Hollenbach et al. (2005) determined the gas mass to 
be < 0.1 Mj between 10-40 AU for a constant surface 
density. 

Later FEPS modeling by Hillenbrand et al. (2008) used 
multi-temperature blackbodies to fit the HD105 disk, and 
found an inner radius of 36.8 AU. This is a lower limit 
to the inner radius since the grains were assumed to be 
large blackbody grains. Smaller grains could reach the 
same temperatures farther from the star. 

The smaller error bars of the Herschel PACS measure- 
ments, and the addition of the 100 /im data point, better 
map the characteristic turnover point in the SED and 



lessens the severity of the model degeneracy. Our results 
are consistent with those of Meyer et al. (2004), Hollen- 
bach et al. (2005), and Hillenbrand et al. (2008). 

The lack of mid-IR excess in this disk indicates an ab- 
sence of both small grains and grains within a large inner 
cavity. The dust begins at about the distance of the So- 
lar System's Kuiper belt (52 AU using our best fitting- 
model) and extends to some distance beyond (see Section 
6 for an estimate of the outer radius of HD105). The lack 
of detectable dust interior to this region may be due to 
one or more planets orbiting inside this dust ring (Moro- 
Marti'n & Malhotra 2005). Apai et al. (2008) searched 
for massive planets in the HD105 disk using VLT/NACO 
angular differential imaging and found none. However, 
the survey was only sensitive to planets with masses > 6 
Mj at distances > 15 AU. 

The minimum grain size for this disk (8.9 /im) is much 
larger than the expected blowout size of 0.5 /im. This 
suggests either the calculation of the blowout size is in- 
accurate, or the small grains are efficiently removed by 
some other mechanism. The calculation of the blowout 
size is just an estimate with several assumptions built 
in, such as grain density and composition, and a change 
in the assumptions could lead to a drastically different 
blowout size. For instance, if the grains were porous, 
they would be less massive at a certain size, and the 
larger grains could be more easily blown out of the sys- 
tem. 

4.2.2. HD202917 

HD202917 is a G5V star 46 pc away. The disk was ob- 
served with Spitzer by Bryden et al. (2006) and analyzed 
as part of the FEPS survey by Hillenbrand et al. (2008) 
and Carpenter et al. (2008). Smith et al. (2006) fit the 
disk with a blackbody grain model to get a lower limit 
on the disk inner radius. They found the inner radius 
to be > 7.4 AU. Hillenbrand et al. (2008) fit a multi- 
temperature blackbody model to the disk that gave an 
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TABLE 5 
Best Fitting Model Parameters 



Target 


n a 






K a 


r a 
' mm 


'F'max a 


g a 


Q b 


Mass (M ( 


b) c 


HD105 


8.9+^/ iiva 


0.5 fim 


1000 fim 


o+U.3 
o.o_ 4 


52+fAU 


120 AC/ 


-1.5 


0.15 


4.3 x 10" 


-4 


HD202917 


< 2.8/j,m d 


0.3 ^m 


1000 fim 


-3.4 ±0.1 


46±Sj AU 


SO AC/ 


-1.5 


0.70 


3.4 x 10" 


-4 


HD3003 


Q K + 0.5 

3-5I j/im 


3.4 fim 


1000 fim 


_4 4+0.1 


7.8l°;^ AU 


120 AC/ 


-1.5 


8.0 x lO" 11 


7.0 x 10" 


-6 


HD12039 


0-4 fim 


0.4 fim 


1000 fim 


-5.5 


14 ± 3 AU 




-1.5 


0.88 


1.7 x 10" 


-4 


HD1466 


0.5 [im 


0.5 fim 


1000 fim 


-3.5 


7.8±?-| AU 


120 AC/ 


-1.5 


0.30 


6.2 x 10" 


-5 



Note. — Model parameters in italics are fixed. 
a Model parameters 
b Model goodness-of-fit 

c Lower limit on dust mass is calculated from the model parameters 
d Unconstrained parameter, limits given are 3cr confidence 



inner radius of 2.5 AU, also a lower limit. HD 202917 
was also resolved in scattered light by HST/ACS, giving 
a disk outer radius of R « 80 AU (Krist 2007; Mustill & 
Wyatt 2009). For this reason, we fixed r max = 80 AU 
rather than the usual value of 120 AU. 

The fit is not well constrained in a m i n , but the fit with 
the lowest \ 2 value has a m i n = 0.3 /jm, equal to the 
blowout size. We found the best fitting inner radius to 
be 46 AU. This result is larger than the previous lower 
limits due to the presence of small grains in the model. 
We find this disk to be consistent with a belt of material 
between 45 and 80 AU. 

4.2.3. HD3003 

HD3003 is an A0V star with a warm disk first detected 
by IRAS (Oudmaijer et al. 1992). Smith et al. (2006) 
also observed the HD3003 disk with Spitzer MIPS at 24 
and 70 fim. Smith et al. (2006) modeled the disk with 
blackbody grains at a single radius and found a dust tem- 
perature of 230 K at a radius of 6.7 AU. Smith & Wy- 
att (2010) added unresolved ground based mid-IR pho- 
tometry from TIMMI2, VISIR, Michelle and TReCS and 
found a blackbody temperature of 265 K with a radius 
of 4 AU. Smith & Wyatt (2010) also make the point that 
HD3003 is a binary, and a disk of this temperature must 
be circumstellar not circumbinary in order to be stable. 

With the new PACS data, we found a lower blackbody 
temperature of 208 K. This is still highest temperature of 
all the disks in the sample, implying the grains are cither 
very close to the star, or smaller than 1 fim in size. But 
HD3003 is an AO star, and has a blowout size of 3.4 /im. 
Therefore, the temperature is likely due to the distance 
from the star. Our model of the HD3003 disk gives an in- 
ner radius of 7.8 AU and a minimum grain size consistent 
with the estimated blowout size. The goodness-of-fit is 
very small for HD3003, with Q = 8 x 10" 11 , despite re- 
producing the PACS data quite well. This small value is 
driven mostly by the small error bars on the Spitzer IRS 
spectrum. However, this goodness-of-fit is indeed larger 
than the value found for the blackbody models. 

The grain size distribution, however, departs from the 
expected shape. The best fit grain size distribution 
power-law index value of k = —4.4 is much steeper than 
the k = —3.5 value expected for a steady state colli- 
sional cascade. HD3003 is the only binary system with a 
confirmed disk in the sample. The apparent binary sep- 
aration is ~ 0.1" (Mason et al. 2001). If this projected 
separation were a true binary separation, the companion 
would be about 4.6 AU from HD3003. Smith & Wyatt 
(2010) suggest the binary must have a semi-major axis 



> 14.4 AU for the disk to be stable. An unstable system 
could provide one explanation for the large departure of 
k from the steady state value. Other possibilities are 
explored in Section 6.3. 

4.2.4. HD12039 

HD12039 is a G3/5V star 42 pc away with a debris disk 
first detected using data from the FEPS Legacy Program 
(Hines et al. 2006). The FEPS team analyzed Spitzer 
observations with IRAC, IRS, and MIPS (24, 70, and 
160 fim). The disk was not detected with Spitzer at 70 
and 160 fim. A blackbody fit gave a characteristic disk 
temperature of 110 K and a lower limit on the radius of 
6 AU. Hines et al. (2006) also fit the disk with the mod- 
els of Wolf & Hillenbrand (2003) that uses astronomical 
silicates, with flat surface density profile (q = 0) and a 
Dohnanyi (1969) distribution of grains from 0.4-1000 fim. 
The best fitting model yielded an inner radius of 28 AU. 

The best-fitting model in Hines et al. (2006) was very 
dependent on the MIPS 70 fim upper limit. With the 
PACS 70 ^m detection, the disk radius is better con- 
strained. We fit the disk with a surface density profile of 
q = —1.5, a minimum grain size at the blowout limit, and 
a Dohnanyi (1969) grain size distribution, similar to the 
modeling done by Hines et al. (2006). But, in contrast 
to Hines et al. (2006), we found an inner radius of 14 ± 3 
AU. This difference is due mainly to our measured disk 
flux at 70 fim that Hines et al. (2006) did not have. 

4.2.5. HDI466 

HD1466 is an F9V star 41 pc away with an excess 
detected by Smith et al. (2006) at both 24 and 70 ^m 
with Spitzer MIPS. Smith et al. (2006) found a minimum 
radius of 7.2 AU. Our physical disk model fit gave an 
inner radius of 7.8 AU, consistent with previous results. 

4.2.6. HD30051 

The HD30051 disk was first discovered in this survey. 
Unfortunately, the disk system was never observed with 
Spitzer, and very little mid-IR data are available to con- 
strain the disk parameters. With only two data points 
showing IR excess, we were unable to model the disk 
in any more detail than a blackbody model. But the 
best-fitting blackbody temperature of 58 K indicates the 
main component of the disk is far from the star. Keeping 
the assumption of large pure blackbody grains, the disk 
would be centered around ~ 45 AU. This is a lower limit 
on the disk radius because smaller silicate grains would 
have a temperature of 58 K tens of AU farther out. 
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Fig. 6. — (a): Azimuthally averaged radial brightness profiles for HDf05 and the reference star Arcturus at 70 /an. The radial profiles 
are azimuthally averaged over annuli with one arcsecond widths, (b): Averaged brightness profiles from contours for HD105 and Arcturus 
at 70 fim. HD105 is extended beyond the PSF in both profiles. 

contours rather than annuli. Figure 7 shows the PACS 
70 /im images of HD105 and the reference star Arcturus 
scaled to the same peak brightness with brightness con- 
tours overlaid. The fluxes in the regions between the 
contours were averaged and the uncertainty measured 
from the standard deviation. The contours used for Fig- 
ure 6 (b) are 15, 30, 50, 70, 85, and 100% of the peak 
brightness for HD105 and Arcturus. For Arcturus, which 
is brighter than HD105, the 1, 2, 5, and 10% contours 
are also used. The radius adopted for each mean flux 
value was the mean distance of the region between the 
contours. This method is unable to map the brightness 
profile of HD105 far from the star because it quickly falls 
below the background. But it it clear from Figure 6 (b) 
that the disk is extended with a FWHM of 7.2 ± 0.5" 
beyond the PSF that has a FWHM of 6.1 ± 0.4". 

The PACS PSF size is given by the Herschel PACS ICC 
as a two-dimensional Gaussian with FWHM of 5.5 x 5.8". 
The FWHM reported here for HD105 was not calculated 
in the same way. Therefore, the FWHM values should 
not be directly compared to the PACS PSF FWHM when 
determining the physical size of the disk. The FWHM 
values are given only to show that the disk is indeed 
significantly extended beyond the PSF. In Section 5.2, 
we determine the physical size of the disk without using 
the FWHM. 



Fig. 7. — Images of HD105 (left) and a PSF reference star Arc- 
turus (right) at 70 fj,m. Brightness contours are overlaid at 15, 30, 
50, 70, 85 and 95% of the peak brightness. 

5. RESOLVING THE HD105 DEBRIS DISK 



5.1. Radial Profile 

One disk in our sample, HD105, is marginally resolved 
at 70 /im in the Herschel PACS images. HD105 is the 
brightest disk in our sample. It may not be the most 
extended, but for the other disks an extended structure 
would fall below the background noise. We determined 
that the HD105 disk was resolved in two ways. The first 
was to compare the azimuthally averaged radial profile 
of the disk to a reference star, Arcturus. Figure 6 (a) 
shows the azimuthally averaged radial profiles of the two 
images. These profiles were calculated by measuring the 
mean flux within annuli of one arcsecond width around 
the star. Scan maps with a smaller pixel scale of 1" 
pixel -1 were created for this purpose. The uncertainties 
in the measurements were determined by calculating the 
standard deviation in the flux within each annulus. This 
method leads to an overestimation of the errors because 
the PACS PSF changes radially across each annulus due 
to its tri-lobed shape. The azimuthally averaged radial 
profile of HD105 shown in Figure 6 (a) is extended be- 
yond the PSF with a FWHM of 6.1 ± 0.8" compared to 
the PSF FWHM of 5.3 ± 0.5". The error in the FWHM 
is determined from the uncertainties in the brightness 
by calculating the distance away from the FWHM where 
the brightness profile plus and minus the errors would 
be equal to half the maximum brightness. This is not 
a strong detection of extended structure given the large 
(and overestimated) errors. Therefore, we also used a 
second method to determine if the disk is resolved. 

The second method tries to avoid problems with the 
PACS PSF tri-lobe shape by averaging over brightness 



5.2. Determining the Outer Radius 

Without data at wavelengths longer than 160 /jm, the 
outer radius cannot be determined from SED fitting 
alone. However, resolved imaging provides geometrical 
information. The marginally resolved image of HD105 
constrains the outer radius of the disk. 

To determine the outer radius of the disk, we pro- 
duced a synthetic image of the disk using the parameters 
from the best fitting SED model. We then convolved 
our model with the PACS PSF 14 and rotated it to the 
position angle of the observation for direct comparison. 
We varied the outer radius and determined the x 2 value 
between the model and the images of HD105. Once the 
outer radius of the model is greater than a certain value, 
the surface brightness at that radius is so low that it 
falls below the noise level. Therefore, we can only place 
a lower limit on the outer radius of the disk. We found 
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a disk outer radius of > 280 AU. 

There are, of course, several assumptions that go into 
this model. The first is that the ring is circular and 
face-on. The second is that the dust surface density dis- 
tribution follows an r -1 - 5 profile. This is expected to 
be the case for collisionally dominated disks. Several 
resolved images of disks have shown more complicated 
structure such as rings with sharper density profiles, 
clumps, warps, and other asymmetries. These structures 
would have less of an effect on the SED, which is fairly in- 
sensitive to the density profile, but would strongly affect 
the analysis of the Herschel images. The outer radius 
derived here is the first such constraint put on the disk, 
but an image with higher spatial resolution is needed to 
better determine the spatial extent. 

6. DISCUSSION 

We detected IR-excesses in ~ 1/3 of the stars in our 
sample. The Tucana-Horologium disks we detect all have 
some similarities. They are optically thin debris disks 
with relatively low fractional luminosities compared with 
younger protoplanetary disks, which have typical frac- 
tional luminosities of Lir/L„ ~ 0.1 (Cieza et al. 2012). 
As of yet, none of the Tucana-Horologium systems have 
shown any significant amount of gas. However, these 
disks do show a remarkable amount of variety for sys- 
tems of the same age (~ 30 Myrs). The single temper- 
ature blackbody fits give a range in temperature of 50- 
208 K. Although the hottest disk is around an AO star 
(HD3003), the other six disks still display a temperature 
range of 50-146 K with no apparent dependence on spec- 
tral type. Any trends that may exist would be hard to 
see in a such a small sample. Our further modeling also 
shows a variety in other disk properties. Three of these 
properties, minimum grain size, the grain size distribu- 
tion power-law, and inner holes, will be discussed in the 
following sections. 

6.1. Minimum Grain Size 

The minimum grain size in debris disks is expected 
to be approximately equal to the blowout size due to 
radiation pressure. Any grains smaller than this limit 
would be ejected from the system on a timescale of a 
few thousand years. But there are a few problems with 
using the blowout size as a limit. First, there are many 
assumptions that go into the blowout size calculation. 
The grains are assumed to be spherical, have a constant 
density, and start out on circular orbits. Additionally, 
the grain composition must be assumed. As circumstellar 
grain composition is hard to determine due to the lack of 
mid-IR solid state features from most debris disks, (e.g. 
Jura et al. 2004), ISM grain composition is commonly 
used. But densities can differ greatly between silicates, 
graphite and ices, all of which are expected to be present 
in debris disks. Porosity may also play a role in changing 
this grain size limit (Lcbrcton et al. 2011). For these 
reasons, it is difficult to interpret the minimum grain 
size results based on the blowout size. 

Only two disks in our sample had well constrained min- 
imum grain sizes. For HD105, the minimum grain size 
found was more than an order or magnitude greater than 
the blowout size. Here, either the calculated blowout size 

14 http://pacs.ster.kuleuven.ac.be/pubtool/PSF/ 



is inaccurate, pointing to a different grain composition or 
porosity, or larger grains are efficiently removed by some 
other mechanism. Thebault & Wu (2008) show that dy- 
namically cold disks have fewer small grains. The smaller 
velocity of the collisions decreases the production rate 
of the small grains while the destruction rate increases. 
This produces a disk dominated by larger grains with or- 
bits mostly confined to the planetesimal ring, very simi- 
lar to what is observed around HD105. For HD3003, the 
minimum grain size is similar to the blowout size. 

6.2. Inner Holes 

The range of temperatures seen in the blackbody fits 
imply inner gaps with radii from 4.5 to 52 AU. The 
largest inner hole is in the HD105 disk. Giant planets 
could be responsible for these holes. However, there are 
also viable mechanisms to explain the holes without plan- 
ets (Kennedy & Wyatt 2010). For instance, dust parent 
bodies may preferentially form at outer locations. This 
can be expected as a result of photoevaporative clearing 
of the inner gas disk and pile-up of dust at its inner edge 
(e.g. Alexander & Armitage 2007) or as a consequence 
of rapid planetesimal formation in spiral arms of a self- 
gravitating disk (e.g. Rice et al. 2006). Resolved imaging 
showing the sharpness of the hole's edge could put more 
constraints on the processes that created the inner hole, 
as was done in the case of the Fomalhaut debris ring 
(Chiang et al. 2009). 

6.3. An Unusual Debris Disk? 

The disk of HD3003 was found to have a grain size dis- 
tribution much steeper than the typical Dohnanyi (1969) 
steady state collisional cascade. A steeper distribution 
implies two possibilities: the disk has an overabundance 
of small grains, or a paucity of large grains. The min- 
imum grain size we find for HD3003 (a m i n = 3.5 /im) 
is large enough that we suspect the explanation to be 
the paucity of large grains. We believe the minimum 
grain size to be correct, as it is consistent with both the 
blowout size and the lack of a 10 /im silicate feature in 
the IRS spectrum which is only seen in disks with small 
(sub-/im) grains (Kessler-Silacci et al. 2006). 

We investigated other explanations of the SED slope. 
The lack of flux at longer wavelengths could mean a 
smaller outer radius than was assumed. The disk could 
be tidally truncated by its binary companion if the other 
star were close enough to the disk. We tried fixing the 
grain size distribution power- law index to k = —3.5 and 
varying r max from 15-120 AU. These models were not 
able to reproduce the observed data. 

The unusual behavior of the HD3003 SED has been 
seen in older debris disks in the DUNES sample (Eiroa 
et al. 2010). Ertel et al. (2012) detected three debris 
disks with Herschel whose slopes are also inconsistent 
with a Dohnanyi (1969) distribution. They give three 
possible explanations for the underabundance of large 
grains. First, there is a departure from the steady state 
collisional cascade conditions. HD3003 is a denser disk, 
so it should be collision dominated not transport domi- 
nated. However, the distribution can still deviate from 
Dohnanyi (1969) by processes such as radiation pressure 
that causes a wavy distribution that is steeper at some 
points (Campo Bagatin et al. 1994; Thebault et al. 2003). 
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The second possibility is that grains of a different compo- 
sition would have a different absorption coefficient, Q a bs, 
which may affect the SED slope. And lastly, a shepard- 
ing planet could result in a spatial separation of the small 
and large grains, leaving the large grains undetectable. 

An idea not proposed in Ertel et al. (2012), but pos- 
sibly relevant here, is enhanced stirring of the planetes- 
imal disk by the companion star. The true separation 
of HD3003's companion star is unknown. If it is close 
enough, it could violently stir the disk. If the orbit of 
the binary is eccentric, then a close passage of the com- 
panion would excite the disk, raising the mean eccen- 
tricity of the disk particles, and thereby increasing their 
relative velocities (Mustill & Wyatt 2009). Numerical 
simulations show that the waviness of the grain size dis- 
tribution depends on the collision velocities. The am- 
plitude and the peak-to-peak wavelength of ripples in 
the grain size distribution increase with higher velocities 
(Thebault et al. 2003; Krivov et al. 2006; Wyatt et al. 
2011). Larger ripples mean the size distribution will be 
steeper for grain sizes above the blowout limit, as is seen 
in HD3003. However, the strange behavior of the SED 
is not yet understood, and will require more data and 
further modeling to determine its cause. 
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TABLE 6 Archive Data Used In SED Fitting. 



HI VI OK 


System 


Wavelength 
W 


Flux Rclcrcncc 
(mjy) 



Hipparcos 

2MASS 

2MASS 

2MASS 

Spitzer/IRAC 

Spitzer/IRAC 

Spitzer/IRAC 

Spitzer/IRAC 

AKARI 

Spitzer/YRS 

Spitzer/MlPS 



0.55 
1.25 
1.65 
2.17 
3.6 
4.5 
5.8 
8.0 
9 
5-37 
24 



3504 ± 32.27 
4139 ± 76.25 
3425 ± 69.41 
2383 ± 43.90 
1023 ± 7.36 

645.4 ± 7.87 

410.5 ±4.27 
230.7 ± 1.52 
223.1 ± 9.59 

28.29 ± 0.25 



H0g et al. 
Cutri et al 
Cutri et al 
Cutri et al 



2000 
2000 
(2003 
(2003 
(2003 



Carpenter et al. (2008) 



Carpenter et al. 
Carpenter et al. 
Carpenter et al. 
Yamamura et al 
Carpenter et al. 
Carpenter ct al. 



2008) 
(2008) 
(2008) 
(2010) 
(2008) 
(2008) 



HmO'2917 


System 


Wavelength 
W 


Flux 

(mjy) 


Rclcrcncc 



Hipparcos 
Hipparcos 
2MASS 
2MASS 
2MASS 
WISE 

Spitzer/IRAC 
Spitzer/IRAC 
WISE 

Spitzer/IRAC 

AKARI 

IRAS 

WISE 

WISE 

Spitzer/IRS 

Spitzer/MIPS 



0.44 
0.55 
1.25 
1.65 
2.17 
3.4 
3.6 
4.5 
4.6 
8.0 
9 
12 
12 
22 
5-37 
24 



619.2 ± 10.84 
1233 ± 14.76 
1770 ± 34.25 
1585 ± 55.47 
1150 ±23.31 

575.86 ± 17.23 

519.2 ± 3.74 
320.8 ± 3.91 
304.36 ± 5.94 

117.3 ± 1.44 
96.55 ± 12.4 
101.5 ± 23.8 
58.93 ± 1.04 
19.82 ± 1.01 

20 ±0.8 



H0g et al 
H0g et al 
Cutri et al 
Cutri et al 
Cutri et al. 
Wright et al 



(2000) 
(2003) 
2003) 
2003) 
(2010) 



Carpenter ct al. (2008) 
Carpenter et al. (2008) 
Wright ct al. (2010) 
Carpenter et al. (2008) 
Yamamura et al. (2010) 
Moshir et al. (1992) 
Wright et al. (2010) 
Wright ct al. (2010) 
Carpenter ct al. (2008) 
Smith ct al. (2006) 



HDM0IW 


Instrument 


Wavelength 
(lim) 


Flux 
(mjy) 


Rclcrcncc 



(2000) 
(2000) 
2003) 
2003) 
2003) 



Hipparcos 

Hipparcos 

2MASS 

2MASS 

2MASS 

AKARI 

IRAS 

AKARI 

Spitzer/IRS 



0.44 
0.55 
1.25 
1.65 
2.17 
9 
12 
18 
5-37 
24 



35400 ± 456.5 
35290 ± 292.6 
15070 ± 513.6 
8870 ±621.4 
6760 ± 124.5 
586.1 ± 9.77 
446.0 ±31.22 
237.7 ± 15.7 



223.9 ±9.0 



Hog et al 
H0g et al 
Cutri et al. 
Cutri et al. 
Cutri et al. 



Yamamura et al. (2010) 
Moshir et al. (1992) 
Yamamura et al. (2010) 
This Work 
Smith et al. (2006) 





Instrument 


Wavclcngti 


1 Flux 


Reference 






(mjy) 





H0g et al. (2000) 
H0g et al. (2000) 
Cutri et al. (2003 
Cutri et al. (2003 
Cutri et al. (2003 
Carpenter et al. 
Carpenter et al. 
Carpenter et al. 
Yamamura et al 
Carpenter et al 
Carpenter ct al 



Hipparcos 
Hipparcos 
2MASS 
2MASS 
2MASS 
Spitzer/IRAC 
Spitzer/IRAC 
Spitzer/IRAC 
AKARI 
Spitzer/IRS 



0.44 
0.55 
1.25 
1.65 
2.17 
3.6 
4.5 
8.0 
9 
7-37 
24 



1103 ± 17.27 
2073 ± 22.91 
2885 ±61.11 
2445 ±85.59 
1718 ±41.15 

747.3 ± 5.38 
470.9 ±5.75 

170.4 ± 1.13 
159.9 ±6.9 

25.65 ±0.23 



2008) 

2008) 

2008) 

(2010) 

2008) 

2008) 



H 1)1466 


Instrument 


Wavelength 
W 


Flux 
(mjy) 


Rclcrcncc 


Hipparcos 
Hipparcos 
2MASS 


0.44 
0.55 
1.25 


2248 ±31.06 
3686 ± 37.34 
4147 ±68.75 


H0g et al. (2000) 
H0g et al. (2000) 
Cutri ct al. (2003) 
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TABLE 6 — continued from previous page 



Instrument 


Wavelength 


Flux 


Hclcrcncc 








(fim) 


(mJy) 








2MASS 


1.65 


3244 ± 107.6 


Cutri et al. 


(2003) 


2MASS 


2.17 


2314 ±36.23 


Cutri et al. 


(2003) 


AKARI 


9 


207.1 ± 20.8 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


141 ± 16.92 


Moshir ct al. (1992) 


Spitzer/IRS 


5-37 




This Work 






Spitzer/MIPS 


24 


32.90 ± 1.3 


Smith ct al 


. (2006) 


H 1 J.iOOfil 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








W 


(mJy) 








Hipparcos 


0.44 


3660 ± 50.56 


H0g et al. ( 


2000) 


Hipparcos 


0.55 


5163 ±47.55 


H0g et al. ( 


2000) 


2MASS 


1.25 


4867 ± 89.66 


Cutri et al. 


(2003) 


2MASS 


1.65 


3829 ± 119.9 


Cutri et al. 


(2003) 


2MASS 


2.17 


2599 ± 52.66 


Cutri et al. 


(2003) 


WISE 


3.4 


1329 ±61.35 


Wright et al. ( 


2010) 


WISE 


4.6 


717.4 ± 14.69 


Wright ct al. ( 


'2010) 


AKARI 


9 


237.4 ±21.9 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


152 ± 15.2 


Moshir et al. ( 


1992) 


WISE 


12 


126.2 ± 2.11 


Wright et al. ( 


'2010) 


WISE 


22 


34.67 ± 1.37 


Wright ct al. ( 


'2010) 


H 1)1 6978 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








(Mm) 


(mJy) 








Hipparcos 


0.44 


95730 ± 1234 


H0g et al. ( 


2000) 


Hipparcos 


0.55 


87120 ± 722.1 


H0g et al. ( 


2000) 


2MASS 


1.25 


26620 ± 7349 


Cutri et al. 


(2003) 


2MASS 


1.65 


17260 ± 4337 


Cutri et al. 


(2003) 


2MASS 


2.17 


13250 ± 439.5 


Cutri et al. 


(2003) 


AKARI 


9 


1068 ±6.11 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


743.0 ± 66.87 


Moshir ct al. (1992) 


AKARI 


18 


257.7 ± 12.1 


Yamamura 


ct 


al. (2010) 


Spitzer/MIPS 


24 


124.0 ±4.96 


Rebull et a 


I- ( 


2008) 


H I MXX4 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








W 


(mJy) 








Hipparcos 


0.44 


79260 ± 1022 


H0g ct al. (2000) 


Hipparcos 


0.55 


72860 ± 604 


H0g et al. (2000) 


2MASS 


1.25 


21720 ± 5128 


Cutri et al. 


(2003) 


2MASS 


1.65 


13790 ± 965.9 


Cutri et al. 


(2003) 


2MASS 


2.17 


10750 ± 356.6 


Cutri et al. 


(2003) 


AKARI 


9 


895.2 ± 31.0 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


1300 ± 78 


Moshir ct al. (1992) 


AKARI 


18 


188.4 ± 14.8 


Yamamura 


ct 


al. (2010) 


Spitzer/IRS 


5-37 




This work 




Spitzer/MIPS 


24 


101.7 ±8.6 


Smith ct al 


. (2006) 


H I >224392 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








(Mm) 


(mJy) 








Hipparcos 


0.44 


37070 ± 478 


H0g et al. ( 


2000) 


Hipparcos 


0.55 


37710 ±312.6 


H0g et al. ( 


2000) 


2MASS 


1.25 


17320 ± 590.3 


Cutri et al. 


(2003) 


2MASS 


1.65 


10730 ± 306.5 


Cutri et al. 


(2003) 


2MASS 


2.17 


7840 ±151.7 


Cutri et al. 


(2003) 


AKARI 


9 


660.5 ± 18.2 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


477 ±28.62 


Moshir ct al. (1992) 


AKARI 


18 


113 ± 15.59 


Yamamura 


ct 


al. (2010) 


H I )28X5 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








(Mm) 


(mJy) 








Hipparcos 


0.44 


54530 ±703.1 


H0g et al. ( 


2000) 


Hipparcos 


0.55 


59330 ±491.8 


H0g et al. ( 


2000) 


2MASS 


1.25 


29820 ± 8003 


Cutri et al. 


(2003) 


2MASS 


1.65 


21340 ± 4192 


Cutri et al. 


(2003) 


2MASS 


2.17 


15160 ± 502.8 


Cutri et al. 


(2003) 


AKARI 


9 


1291 ± 18.6 


Yamamura 


ct 


al. (2010) 


AKARI 


18 


340.6 ± 25.6 


Yamamura 


ct 


al. (2010) 


Spitzer/MIPS 


24 


156.1 ±4.5 


Smith ct al 


. (2006) 


H I )53X42 


Instrument 


Wavelength 


Flux 


Rclcrcncc 








(Mm) 


(mJy) 








Hipparcos 


0.44 


2529 ± 37.28 


H0g et al. ( 


2000) 


Hipparcos 


0.55 


3716 ± 37.65 


H0g et al. ( 


2000) 


2MASS 


1.25 


3813 ± 101.9 


Cutri et al. 


(2003) 


2MASS 


1.65 


2831 ± 80.83 


Cutri et al. 


(2003) 


2MASS 


2.17 


2015 ± 38.98 


Cutri et al. 


(2003) 


WISE 


3.4 


894.4 ±31.9 


Wright et al. (2010) 


WISE 


4.6 


525.0 ± 10.7 


Wright ct al. (2010) 


AKARI 


9 


179.4 ±8.49 


Yamamura 


ct 


al. (2010) 


IRAS 


12 


140.0 ± 10.0 


Moshir et al. ( 


1992) 


WISE 


12 


100.1 ± 16.7 


Wright et al. ( 


2010) 


WISE 


22 


40.4 ±1.29 


Wright et al. ( 


'2010) 


Spitzer/IRS 


7-37 




Moor et al. 


(2009) 
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Instrument 


Wavelength 
W 


Flux 
(mjy) 


Hclcrcncc 


Spitzer/MIPS 


24 


31.0 ±1.3 


Moor et al. (2009) 


H 1 )44627 


Instrument 


Wavelength 


b'lux 
(mjy) 


Hclcrcncc 



Hipparcos 

2MASS 

2MASS 

2MASS 

WISE 

WISE 

AKARI 

IRAS 

WISE 

WISE 



Instrument 



0.55 
1.25 
1.65 
2.17 
3.4 
4.6 
9 
12 
12 
22 



Wavelength 
(Mm) 



,.21 

767.2 ± 12.72 
1486 ± 32.86 
1497 ± 28.95 
1075 ± 23.77 
533.0 ±14.93 
295.5 ± 5.77 
181 ±5.09 
101 ± 16.16 
58.39 ± 1.03 
15.06 ± 0.81 

3D m 




H0g et al. 
Cutri et al. 
Cutri et al. 
Cutri et al. 
Wright ct al 
Wright et al 
Yamamura et al. (2010) 
Moshir et al. (1992) 
Wright et al. (2010) 
Wright ct al. (2010) 



JmJy]_ 



Hclcrcncc 



H0g et al. (2000) 
H0g et al. (2000) 
Cutri et al. (2003) 
Cutri et al. (2003) 
Cutri et al. (2003) 
Wright et al. (2010) 
Wright et al. (2010) 
Yamamura et al. (2010) 
Wright ct al. (2010) 
Wright ct al. (2010) 



Hipparcos 

Hipparcos 

2MASS 

2MASS 

2MASS 

WISE 

WISE 

AKARI 

WISE 

WISE 



Instrument 



0.44 
0.55 
1.25 
1.65 
2.17 
3.4 
4.6 
9 
12 
22 



Wavelength 
(Mm) 



110.5 ± 6.01 

304.0 ±8.96 
787.9 ± 16.69 
754.9 ±39.65 
579.6 ± 13.88 
273.9 ± 6.38 

148.1 ± 2.89 
65.96 ± 5.85 
29.78 ± 0.61 
8.36 ± 1.23 



lux 

(mJy) 



Reference 



Hipparcos 

Hipparcos 

2MASS 

2MASS 

2MASS 

AKARI 

IRAS 



Instrument 



0.44 
0.55 
1.25 
1.65 
2.17 
9 
12 



Wavelength 



141.1 ±6.239 
512.0 ±9.904 
1852 ±30.71 
2098 ± 65.72 
1625 ± 26.93 
162.7 ± 11.6 
121.0 ±20.57 
HIP107345 — 
Flux 



H0g et al. (2000) 
H0g et al. (2000) 
Cutri et al. (2003) 
Cutri et al. (2003) 
Cutri et al. (2003) 
Yamamura et al. (2010) 
Moshir ct al. (1992) 



Hclcrcncc 



(Mm) 


(mJy) 


0.44 


16.24 ±4.654 


0.55 


78.94 ±8.231 


1.25 


503.6 ± 11.60 


1.65 


596.3 ± 12.08 


2.17 


472.4 ±11.31 


3.4 


236.5 ± 5.51 


4.6 


135.1 ± 2.64 


12 


27.71 ±0.54 


22 


7.88 ±0.89 



H0g et al. (2000) 
H0g et al. (2000) 
Cutri ct al. ( 
Cutri et al. { 
Cutri et al. I 
Wright ct al 
Wright ct al 
Wright ct al 
Wright ct al 



Hipparcos 

Hipparcos 

2MASS 

2MASS 

2MASS 

WISE 

WISE 

WISE 

WISE 



2003) 
2003) 
2003) 
(2010) 
(2010) 
2010) 
2010) 



H 1 P3K56 


Instrument 


Wavelength 


Flux 


Hclcrcncc 






(fim) 


(mJy) 






Hipparcos 


0.44 


17.48 ±4.012 


H0g et al. (2000) 


Hipparcos 


0.55 


41.43 ±5.780 


H0g et al. (2OOO) 


2MASS 


1.25 


645.8 ± 11.90 


Cutri et al. ( 


'2003) 


2MASS 


1.65 


730.3 ± 16.14 


Cutri et al. 1 


'2003) 


2MASS 


2.17 


595.3 ± 14.81 


Cutri et al. ( 


'2003) 


AKARI 


9 


118.0 ± 14.8 


Yamamura ct al. (2010) 


Spitzer/MIPS 


24 


8.4 ±0.34 


Rcbull et al. 


(2008) 


(;KC80Kfi-4K2 


Instrument 


Wavelength 


Flux 


Reference 






(Mm) 


(mjy) 






Zciss/FOTRAP 


0.44 


15.13 ±4.82 


Torres ct al. 


(2006) 


Zeiss/FOTRAP 


0.55 


54.15 ± 9.77 


Torres et al. 


(2006) 


2MASS 


1.25 


683 ± 14.62 


Cutri et al. 1 


'2003) 


2MASS 


1.65 


809.66 ± 17.33 


Cutri et al. { 


'2003) 


2MASS 


2.17 


666.19 ± 16.77 


Cutri et al. 1 


'2003) 


AKARI 


9 


86.52 ± 14.8 


Yamamura c 


t al. (2010) 


Spitzer/MIPS 


24 


9.0 ± 0.36 


Rcbull et al. 


(2008) 



